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A set of universal relations between various properties of any few-body or many-body system 
consisting of fermions with two spin states and a large but finite scattering length have been derived 
by Shina Tan. We derive generalizations of the Tan relations for a two-channel model for fermions 
near a Feshbach resonance that includes a molecular state whose detuning energy controls the 
scattering length. We use quantum field theory methods, including renormalization and the operator 
product expansion, to derive these relations. They reduce to the Tan relations as the scattering 
Ch ' length is made increasingly large. 
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I. INTRODUCTION 



Many-body systems of fermions have long been of great importance in astrophysics, nuclear physics, and solid state 
physics. The development of trapping and cooling techniques for ultracold atoms has made them important in atomic 
d • physics as well. In this case, the strength of the interaction is governed by the 2-body scattering length which can be 
^ I controlled experimentally, adding a new dimension to the problem [1] . 

' If the scattering length a is much larger than the range of the interactions, the system has universal properties that 
are determined only by the large scattering length. For sufficiently low number density n, the universal properties can 
be calculated using perturbative methods. If n\a\^ is comparable to 1 or larger, the problem is nonperturbative. In 
Q ■ the special case of two equally-populated spin states, systematically improvable calculations are possible using Monte 
Carlo methods [2]. If the populations are not equal, this approach suffers from the fermion sign problem. If there are 
three or more spin states, the problem is complicated by the Efimov effect [3]. The challenging nature of the general 
problem makes exact results very valuable. One case in which exact results are known is the unitary limit a = ±oo, 
where they can be derived by exploiting scale invariance [4] and conformal invariance [5] . 
J — ' Shina Tan has derived a number of universal relations between various properties of an arbitrary system consisting 
P\J of fermions in two spin states with a large scattering length [6-8]. The Tan relations include the coefficient of the 
tail in the momentum distribution [6] , a decomposition of the energy E into terms that are insensitive to short 
\^ , distances [6] , an expression for the local pair density [6] , the rate of change of E from changes in the scattering length 
' [7], the relation between the pressure and the energy density in a homogeneous system [7], the virial theorem for a 
00 ^ system in a harmonic trapping potential [8], and an expression for the inelastic 2-body loss rate [9]. These relations 
' all involve a property of the system that Tan called the integrated contact intensity. For brevity, we will refer to it 
simply as the contact and denote it by C. It can be expressed as the integral over space of the contact density., which 
we will denote by C. The Tan relations hold for any state of the system: few-body or many-body, homogeneous or in a 
trapping potential, superfluid or normal, zero or nonzero temperature. Tan derived his relations within the framework 
of the many-body Schrodinger equation. He used novel methods involving generalized functions to implement the 
Bethe-Peierls boundary condition associated with the large 2-body scattering length. 

In Ref. [10], Braaten and Platter derived the Tan relations using a quantum field theory formulation of the problem. 
They identified the contact density C as the expectation value of a local composite operator constructed from the 
quantum fields. To derive the Tan relations, they used standard methods of renormalization together with the operator 
product expansion (OPE). The OPE was invented independently by Ken Wilson [11] and by Leo Kadanoff [12] in 1969. 
The OPE expresses the product of local operators separated by a short distance as an expansion in local operators 
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with coefficients that are functions of the separation of the operators: 

Oa{R - \r)OB{R +\r) = Y, ClsirWdR). (1) 

c 

The sum is in general an infinite sum over all possible local operators Oc{R)- The functions C^^(r) are called 
Wilson coefficients or short- distance coefficients. Some of the Wilson coefficients can be singular at r = 0. The OPE 
is an asymptotic expansion for small r : for any power N, only a finite number of terms on the right side of the OPE 
have Wilson coefficients that go to zero more slowly than as r ^ [13]. 

Wilson proposed the OPE as an alternative framework to the Lagrangian or Hamiltonian formulations of quantum 
field theory [11]. His motivation was the problem of the strong interactions associated with the nuclear force. Some 
aspects of the strong interactions had been successfully explained using current algebra, which is the algebra satisfied 
by the currents associated with the symmetries of a quantum field theory. For example, if the charge operators 
associated with a symmetry satisfy the Lie algebra [Q", = if'^'"^Q'^, the associated charge density operators p"'{R) 
satisfy the equal-time commutation relations 

[p%R - ir)), p'iR + ir)] = ir''-S^r)p%R). (2) 

This is a special case of the general short-distance expansion in Eq. (1). Wilson observed that in a scale-invariant 
theory, the dependence of the coefficients ^ (r) on the separation r is determined up to a multiplicative constant by 
scale invariance and other symmetries. He used this observation to derive nontrivial results for the strong interactions 
that followed only from the assumption of scale invariance. The quantum field theory that describes the strong inter- 
actions was eventually determined to be quantum chromodynamics (QCD). Wilson's assumption of scale invariance 
in the strong interactions was justified by the discovery that the coupling constant of QCD is asymptotically free 
at short distances [14, 15]. QCD and the electroweak gauge theory that together comprise the Standard Model of 
elementary particle physics are Lagrangian quantum field theories. Thus the OPE has not proved to be essential as 
an alternative formulation of the relativistic quantum field theories that describe elementary particles. However the 
OPE has become a standard tool for developing systematic approximations to the Standard Model [16]. 

Kadanoff used the OPE to understand critical phenomena in condensed matter physics [12]. He showed that the 
critical exponents that describe the scaling behavior of correlations finictions at a critical point can be deduced from 
the knowledge of which Wilson coefficients are singular as r ^ 0. The OPE therefore provides powerful constraints on 
critical behavior in condensed matter systems [17]. A critical point is usually characterized not only by scale invariance 
but also by conformal invariance. The long-distance behavior near the critical point can therefore be described by 
a statistical field theory with conformal symmetry, i.e. a conformal field theory. The OPE is one of the basic tools 
that is used to study conformal field theories. The quantum field theory that describes fermionic atoms with two 
spin states in the unitary limit is a conformal field theory. The classification of local operators based on conformal 
symmetry has been exploited by Nishida and Son [18] and by Mehen [19]. However, until Ref. [10], the OPE had not 
been applied directly to the problem of fermionic atoms with a large scattering length. 

The T-matrix element for S-wave atom-atom scattering can be written as 

Att 1 
^ ' " '^kcoi5o{k)-ik' 

where m is the mass of the atoms, k is the relative wavenumber, and (5o(fc) is the S-wave phase shift. If the interactions 
have a finite range, the phase shift has a low energy expansion that is conveniently expressed in the form 

kcoi6o{k) = -l/a+ ^Tsk^ + . . . . (4) 

This expansion, which is called the effective range expansion, defines the scattering length a and the effective range 
. The Tan relations apply to any system in which the phase shift at the accessible energies can be approximated by 

the first term in the effective range expansion. We will consider the generalizations of these relations in two scattering 

models that are specified by the S-wave phase shift. 
In the Zero-Range Model, the phase shift is given by the leading term in the effective range expansion up to 

arbitrarily high energy: 

fccot(5o(fc) = -1/a. (5) 

This model can be formulated as a renormalizable local quantum field theory. Braaten and Platter derived the Tan 

relations by applying standard renormalization methods and the OPE to the quantum field theory for the Zero- Range 
Model [10]. One advantage of this derivation is that it can be generalized in a straightforward way to any model that 
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can be formulated as a renormalizable local quantum field theory [20]. One such model is the Resonance Model in 
which the S-wave phase shift is given by 

kcotSoik) = -A-K (x — 

\ mu — 

where v, g, and A are real parameters. This model was first formulated as a quantum field theory by Kaplan [21] and 
by Kokkclmans ct al. [22] . It is of practical importance in cold atom physics because it provides a natural description 
of fermions near a Feshbach resonance. An important special case of the Resonance Model is the Effective Range 
Model, which can be obtained by setting A = 0. The equation for the S-wave phase shift in Eq. (6) reduces in this 
limit to the form 

k cot So{k) = -l/a+ ^rgk"^, (7) 

where a is the scattering length and < is the effective range. The effective range must be negative for this phase 
shift to arise from a local quantum field theory that is renormalizable [23]. The Effective Range Model has been 
proposed by Petrov as a model for a narrow Feshbach resonance [24] . 

In this paper, we derive the generalizations of the Tan relations for the Resonance Model and the Effective Range 
Model. We begin in Section II by enumerating the Tan relations. In Section III, we summarize the derivation of 
these relations in Ref. [10] using the quantum field theory for the Zero-Range Model. In Section IV, we describe 
the formulation of the Resonance Model as a renormalizable local quantum field theory. The generalizations of the 
Tan relations in the Resonance Model are derived in Sections V-XI. We use the OPE to identify the contact density 
operator. The generalizations of the Tan relations involve the expectation values of two other local operators besides 
the contact density operator. In limits in which the scattering length becomes arbitrarily large, the generalized Tan 
relations of the Resonance Model must reduce to the original Tan relations. The resulting constraints on matrix 
elements of local operators are derived in Section XII. In Section XIII, we derive the generalizations of the Tan 
relations for the Effective Range Model by exploiting the fact that it is a limiting case of the Resonance Model. Our 
results are summarized in Section XIV. The diagrammatic methods that we use to calculate the Wilson coefficients 
in the OPE are described in Appendix A. 




II. THE TAN RELATIONS 



In a scries of papers [6- 8], Shina Tan derived some universal relations for states consisting of fermions with two 
spin states that interact through a large scattering length a and may be trapped in an external potential V{R). We 
will refer to the fermions as atoms, although this problem also has applications in condensed matter physics and in 
nuclear physics. We label the two spin states by cr = 1, 2. The Hamiltonian can be expressed as the sum of a kinetic 
term T, an interaction term /, and an external potential term V: H = T + I + V . The Tan relations all involve a 
quantity C that we call the contact. There are seven distinct Tan relations: 

1. Tail of the momentum distribution [6]. The momentum distribution pa{k) for atoms with spin a has a 
power-law tail that decreases at large wavenumber k like 1/k'^: 

p,{k) C/k\ (8) 

The coefficient C is the same for both spin states. 

2. Energy relation [6]. The energy E = {H) = (T) -|- (/) + {V) can be decomposed into four terms that are each 
insensitive to distances much smaller than \a\: 

^ ^ t d^k P f C\ C 

The interaction energy (7) has been separated into two terms proportional to C. One of them has been combined 
with the kinetic term to make the integral over k convergent. 

3. Local pair density [6]. If A^pair(-R, s) is the number of pairs of atoms with spins 1 and 2 in a small ball of 
volume Itts"^ centered at the point R, the contact density C at that point can be expressed as 

C{R) = lim(4/s4)iVp^i^(i2, s). (10) 

s— »0 



Note that iVpair(i?, s) scales as s** as s ^ instead of s^, so the number of pairs with small separation s is much 
larger than one might naively expect. 
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4. Adiabatic relation [7] . The rate of change of the energy E = {H) due to a change in the scattering length is 

da Anina'^ ^ ^ 

5. Pressure relation [7, 8]. For the homogeneous system with constant external potential, the pressure V and 
the energy density £ are related by 

where C is the contact density. 

6. Virial theorem [8]. The virial theorem for atoms in a harmonic trapping potential in the unitary limit a = ±oo 
was derived in Ref. [25]. Tan derived the virial theorem for the case of finite a [8]: 

E = 2(V) - — — C. (13) 

7. Inelastic 2-body loss rate [9]. If a pair of atoms has inelastic scattering channels, the scattering length a 
has a negative imaginary part. To leading order in Ima, the rate at which the atoms are depleted by inelastic 
collisions is 

r « ^^^^a (14) 



The adiabatic relation in Eq. (11) can be used as an operational definition of the contact. A simple example in the 
case a > is the shallow diatomic molecule, or dimer, formed by atoms with spins 1 and 2. The energy of a dimer at 
rest is E = —h^/{ma'^). Prom Eq. (11), the contact for the dimer is 

C = Sw/a. (15) 

In Ref. [7], Tan used the adiabatic relation in Eq. (11) to determine the contact density for the homogeneous gas 

consisting of equal populations of atoms with spins 1 and 2 in three Umits: the BCS limit (a ^ 0~), the unitary limit 
(a — * ±oo), and the EEC limit (a 0+). If the total number density of atoms is n, the contact densities in these 
three limits are 

as a — > 0~, (16a) 

as a ^ ±00, (16b) 

as a ^ 0^, (16c) 

where ( in Eq. (16b) is a universal constant whose value was estimated in Ref. [7] to be approximately 1. 
The most promising systems for testing the Tan relations experimentally are cold trapped atoms near a Feshbach 

resonance. Some of the possibilities have been discussed by Tan in Refs. [6, 7]. The momentum distributions Paik) 
can be measured by suddenly turning off the trapping potential and simultaneously using the Feshbach resonance to 
change the scattering length to 0. The cloud of atoms will expand and the initial momentum distribution can be 
determined from its density distributions after the expansion. There should be a sealing region in which Pa{k) has 
the form in Eq. (8) and the contact C is simply the coefficient. As pointed out by Tan in Ref. [6], the energy relation 
in Eq. (9) can be used to determine the sum (T) + (/) of the kinetic and interaction energies from the momentum 
distribution. The potential energy {V) can be determined from the density distributions before the expansion. The 
virial theorem in Eq. (13) is a nontrivial relation between (T) + (/), {V), and C that can be tested experimentally. 
The adiabatic relation in Eq. (11) can also be tested by using the Feshbach resonance to change the scattering length 
and then measuring the changes in (T) + (/), {V), and C. 
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III. ZERO-RANGE MODEL 



In Ref. [10], the Tan relations were rederived using the quantum field theory for the Zero-Range Model. They 
follow straightforwardly from standard renormalization methods and from the operator product expansion (OPE). 
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This analysis revealed that the contact C can be expressed as the integral over space of the expectation value of a 
local operator. In this section, we summarize the derivations in Ref. [10]. 

A quantmn field theory that describe atoms with two spin states must have two quantum fields V'o-l'"), cr — 1,2. 
The Hamiltonian for a local quantum field theory can be expressed as the integral over space of a Hamiltonian density: 
H = J (Pr Ti. If the atoms are in an external potential V{r), the Hamiltonian density is the sum of a kinetic term T, 
an interaction term X, and an external potential term V: 7i = T + X + V. The simplest quantum field theory that 
describes atoms with a large scattering length is the Zero-Range Model, in which the phase shift has the simple form 
given in Eq. (5) up to arbitrarily large momentum. For the Zero-Range Model, the three terms in the Hamiltonian 
density are 

cr 

I = ^4i,U,4^\ (17b) 
V = V{R)Y,€i^.- (17c) 

For simplicity, we have set h = 1. The superscripts (A) on the operators in Eqs. (17a) and (17b) indicate that their 
matrix elements are ultraviolet divergent and an ultraviolet cutoff is required to make them well defined. For the 
ultraviolet cutoff, we impose an upper limit |fe| < A on the momenta of virtual particles. In the limit A ^ oo, the 
Hamiltonian density in Eq. (17) describes atoms with the phase shift given by Eq. (5) if we take the coupling constant 
to be 

The connected Green's function for the scattering of a pair of atoms with spins 1 and 2 has a well-behaved limit as 
A ^ oo: 

A{E) = . , (19) 

— 1/a + y'—mE — le 

where E is the total energy of the two atoms in their center-of-mass frame. The T-matrix element for scattering of a 
pair of atoms with momenta +p and —p is obtained by setting E = p^/m. Comparing with the expression for the 
T-matrix element in Eq. (3), we obtain the phase shift in Eq. (5). 
In Ref. [10], the contact C in the Tan relations was identified as the integral over space of the expectation value of 

the local operator Aq(A)')/'iV'2V'i''/'2^''- Matrix elements of the operator, which is the product of four quantum fields, 
are ultraviolet divergent, as indicated by the superscript (A), but the dependence of the matrix elements on A is 
precisely cancelled by the prefactor Aq(A). Since this operator has ultraviolet finite matrix elements, we suppress the 
dependence on A and denote it by AgV'iV'IV'iV's- The contact can be expressed as^ 

C = Jd^R{X\Xli;l4i;iMR)\X). (20) 

This expression makes explicit the dependence of the contact on the state \X). It also reveals that the contact is an 
extensive quantity. 

For some of the Tan relations, the derivations in Ref. [10] followed in a straightforward way from the renormalization 

of the quantum field theory once one realizes that Aq^/j|7/'|'0iV-'2 is a finite operator. For example, the energy relation 
in Eq. (9) can be derived by separating the interaction term / d'^R{X\l\X) into two terms by using the identity 

The adiabatic relation in Eq. (11) can be derived from the Feynman-Hellman theorem, 

dE/da = jd^R{X\dH/da\X), (22) 



^ We denote the product of local operators Oa{R) and C's(H) at the same point in space by OaOb{R)- 
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together with the identity 



The pressure relation in Eq. (12) can be derived by applying dimensional analysis to the free energy density (or 
thermodynamic potential density) which requires 



9__1 d_ 



T=\t. (24) 



The virial theorem in Eq. (13) can be derived by applying dimensional analysis to the energy E = f (fR{X\'H\X), 
which requires 



' d Id 
doj 2 da 



E = E. (25) 



The expression for the inelastic loss rate F in Eq. (14) can be derived from the adiabatic relation in Eq. (11) by 
identifying — as the imaginary part of the energy E that results from adding a small negative imaginary part to 
the scattering length a to take into account the effects of the inelastic scattering channels. 

The derivations of the tail of the momentum distribution in Eq. (8) and the expression for the contact density C 
in terms of the local pair density in Eq. (10) require the OPE. The number density operator for atoms with spin a is 
■tp^tjjaiR)- The momentum distribution for atoms with spin a in the state \X) can be expressed as 

p^ik) = {X\i>lik)Mk)\X), (26) 

where tpa{k) is the Fourier transform of the quantum field il>a{r). The expression for the momentum distribution in 
terms of the quantum field is 

Pa{k) = jd^'R J(fre^''-'-{X\ijl{R-^r)MRHr)\X). (27) 

The OPE for the bilocal operator inside the expectation value has the form 

iji{R-^r)MR+ ir) =Y,Ca,n{r)On{R), (28) 

n 

where the sum is over all possible local operators. The local operators can be expressed as products of any number of 

quantum fields ipa{R) and an equal number of quantum fields ^jj'l{R), with any number of gradients applied to those 
fields. If the OPE in Eq. (28) is inserted into Eq. (27), the momentum distribution reduces to 

Paik) = ^Q d\ e^'=-'-a,„(r)) j d^R {X\Or^{R)\X). (29) 

Since the OPE is an asymptotic expansion for small r, Eq. (29) is an asymptotic expansion for large fc. If a Wilson 
coefficient Ccr,n(r) can be expanded as a power series in the vector r, the corresponding integral in Eq. (29) can be 
expressed in terms of the Dirac delta function in k and derivatives of the Dirac delta function. A power-law tail in 
Pa{k) at large k can arise only from a term whose Wilson coefficient is not an analytic function of the vector r at 
r = 0. 

As shown in Ref. [10], the OPE in Eq. (28) can be expressed as 

i>l{R- ir)V.(H+ ir) = i^lMR) + \r ■ [V4VV.(i2) - VViV.(i2)] 

-^A^VlV'IV'iV'2(i?) + ..., (30) 

where we have written explicitly all terms whose Wilson coefficients go to zero more slowly than as r — > 0. The 
first two terms on the right side of the OPE in Eq. (30) can be obtained by multiplying the Taylor expansions of the 
two operators. The third term arises from quantum fluctuations involving pairs of atoms with small separations. Its 
Wilson coefficient is proportional to r = |r|, which is not an analj^tic function of the vector r at r = 0. Its Fourier 
transform at nonzero values of k is given by 



/ 



rf3rre*'"^ = -|j. (31) 
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This can be derived by differentiating the Fourier transform of 1/r. The corresponding term in Eq. (29) gives a 
power-law tail in the momentum distribution: 

p^k)^ Jd^R{X\Xl44ij,MR)\X)/k\ (32) 

Comparing with Eq. (8), we obtain the expression for the contact C in Eq. (20). 

The expression for C in terms of the local pair density in Eq. (10) follows from the OPE for ■^|i/»i(ii + r) and 
tpltP2{R + r'). In Ref. [10], it was shown that the OPE includes the operator Aq'^|'^|'^iV'2 with a Wilson coefficient 
proportional to \r' — r\~^: 

4MR + r)4MR + r') = ^Q^2y _ ^^2 >^l44i^iMR) + ■■■■ (33) 

Only the most singular term as r' ^ r is shown explicitly on the right side. Integrating the left side over r' and r 
inside a small ball of radius s, wc obtain an operator that counts the number of pairs of atoms with spins 1 and 2 
inside that ball. If the right side of Eq. (33) is integrated over the same region and multiplied by An/s, the only term 
that survives in the limit s — > is the one that is shown explicitly. Taking the expectation value of both sides, we 
obtain the expression for the contact density C in Eq. (10). 



IV. RESONANCE MODEL 



A practical method for controlling the scattering length of atoms is by exploiting Feshbach resonances. A Feshbach 

resonance occurs at a value Bq of the magnetic field for which a molecular state in a closed channel with a higher 
two-atom threshold is in resonance with a pair of atoms at threshold in the channel of interest. Near the Feshbach 
resonance, the dependence of the scattering length a on the magnetic field B can be approximated by 



a{B) = abg (1 - -g— ^ J , (34) 

where abg is the background scattering length far from the resonance and -h A is the position of a zero of the 
scattering length. 

A natural description for atoms near a Feshbach resonance is provided by the Resonance Model, in which the 
molecule responsible for the resonance is treated as a point particle. The Resonance Model has three parameters ly, 
g, and A that can be defined by the S-wave phase shift given in Eq. (6). The scattering length and the effective range 
are 

a = f A - ^) , (35a) 



4n \ mv 

= -87r(x-^] (35b) 

The effective range is negative definite. The scattering length in Eq. (35a) can be made arbitrarily large by tuning 
ly to near 0. In this case, the limiting value of the effective range is —Sn/g'^. The standard expression for the 

scattering length near a Feshbach resonance in Eq. (34) can be reproduced by taking u to be linear in the magnetic 
field B while g^ and A are constants: 

p = -n{B-Bo), (36a) 
g^ = —iTTmahg fi A, (36b) 
A = 47rabg. (36c) 

The parameter fj, in Eqs. (36) is the difference between the magnetic moment of the molecule and twice the magnetic 
moment of an isolated atom. Away from the Feshbach resonance, the parameter v can be interpreted as the detuning 
energy of the molecular state relative to the two-atom threshold. Near the Feshbach resonance, the strong coupling 
between the molecular state and two-atom scattering states makes the interpretation of i/ more complicated. 

The Resonance Model was first formulated as a local quantum field theory by Kaplan [21] and by Kokkelmans et 
al. [22]. In the sector that consists of states containing two atoms or one diatomic molecule, the model can be solved 
analytically. The solution reveals that the model is renormalizable, at least in this sector. Kaplan used dimensional 
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FIG. 1: Integral equation for the scattering of a pair of atoms in the Resonance Model. Single lines denote atom propagators and 
double lines denote molecule propagators. The blob represents the amplitude iA{E). The solid dots represent the interaction 
vertices —iXo/m and —igo/m. 



regularization to remove ultraviolet divergences, thus avoiding the need for explicit renormalization of the parameters 
[21]. Kokkclmans ct al. [22] derived the rcnormahzations of the parameters that arc required to make the obscrvablcs 
independent of the ultraviolet momentum cutoff. They used numerical solutions of mean-field integral equations for 
this model to demonstrate the feasibility for achieving superfluidity in ultracold gas of fermionic atoms. 

In the quantum field theory formulation of the Resonance Model for fermionic atoms with two spin states, there 
are three quantum fields: fermionic fields f/'i and ■i/'2 that annihilate atoms and a bosonic field </> that annihilates a 
diatomic molecule. The Hamiltonian density for the Resonance Model is the sum of a kinetic term T, an interaction 
term I, and an external potential term V: W = T + J + V. The individual terms are 



T 

1 
V 



™ 2m 

F(il)(ViVi+V'IV'2 + 2 



2m 



4m 
Ao 
m 



(37a) 
(37b) 
(37c) 



To avoid ultraviolet divergences, an ultraviolet cutoff A must be imposed on the momenta of virtual particles. The 

superscripts (A) on some of the operators in Eqs. (37a) and (37b) indicate that they have matrix elements that arc 
ultraviolet divergent. For the model to have a nontrivial finite limit as A — > cx), the parameters vq^ dOi and Aq in the 
interaction term in Eq. (37b) must depend on A. We will refer to them as the bare parameters. The parameters u, g, 
and A in the phase shift in Eq. (6) do not depend on A. We will refer to them as the renormalized parameters. The 
values for the bare parameters that are required to reproduce the phase shift in Eq. (6) are 



Ao(A) = Z(A)A, 



where the renormalization constant Z is 



30 (A) : 

1.0(A) : 

Z{A) 



ZiA)g, 



AAy' 

27r2 ) 



There are two simple combinations of the parameters that are renormalization invariants: 



5 
A 

mX 



Ac' 

= 2^0- 



mXo 



(38a) 
(38b) 

(38c) 



(39) 

(40a) 
(40b) 
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(a) (6) 



FIG. 2: Diagrammatic representations of (a) the bilocal operator ipl{R— ^r)'ipa(R+ ^r) and (b) local operators that annihilate 
one atom and create one atom, such as ip^ipaiR)- 



The amplitude for the scattering of a pair of atoms can be calculated by solving the Lippmann-Schwinger integral 
equation, which is represented diagrammatically in Fig. 1. After evaluating the loop integral using Eq. (Al4a), the 
integral equation reduces to an algebraic equation: 



A{E) = --{Xo + 
m 



E — vo + ie 



A 



le 



(41) 



where A is an ultraviolet momentum cutoff. The solution A{E) depends on the total energy E of the pair of atoms in 
the centcr-of-mass frame and not separately on their momenta. After substituting Eqs. (38) for the bare parameters, 
the amplitude becomes independent of the ultraviolet cutoff: 



A{E) 



1 

m 



X + 



E-v 



1 

4-77 



^/—mE — ie 



(42) 



The T-matrix element for scattering of a pair of atoms with momenta +p and —p is obtained by setting E = p'^/m. 
Comparing with the expression for the T-matrix element in Eq. (3), we obtain the phase shift in Eq. (6). 

Since the Resonance Model can be formulated as a quantum field theory, the methods used to derive the Tan 
relations in Ref. [10] can be applied equally well to this model. Because the Resonance Model is a local quantum 
field theory, we can use the OPE to identify the contact density operator. Because it is renormalizable, we can use 
standard renormalization methods to derive the analogs of some of the other Tan relations. In the following seven 
sections, we deduce the generalized Tan relations for the Resonance Model. 



V. TAIL OF THE MOMENTUM DISTRIBUTION 



The Tan relation in Eq. (8) states that the momentum distributions Pa{k) have power-law tails proportional to 
with the coefficient for both spin states given by the contact C. We will use the operator product expansion 
(OPE) to show that the momentum distributions in the Resonance Model also have power-law tails proportional to 
l/fc** and we will express the contact C as the integral over space of the expectation value of a local operator. 

Our starting point is the expression for the momentum distribution Pa{k) in Eq. (27), which involves the expectation 
value of V'J-(-R— |»')V'<t(-R+ 1^)- Inserting the OPE in Eq. (28) for this bilocal operator, we obtain the asymptotic 
expansion of Pcr{k) for large k in Eq. (29). A power- law tail in Pa-{k) can only come from terms in the OPE whose 
Wilson coeflScients Ca^nif) are not analytic functions of r at r = 0. Such terms can arise from short-distance quantum 
fluctuations involving pairs of atoms. 

The Wilson coefficients on the right side of the general OPE in Eq. (1) can be determined by a procedure called 
matching. First, matrix elements of the bilocal operator Oa{R— ■^r)OB{R + ^r) between appropriate states are 
calculated and expanded in powers of r = |r|. Second, matrix elements of the local operators Oc{R) between the 
same states are calculated. Finally, the coefficients g{r) are determined by demanding that the expansions in 
powers of r on both sides of the equation match. The coefficient of Oc{R) can be determined by matching matrix 
elements between any states \X) and \X') for which {X'\OciR)\X) is nonzero. We can exploit the fact that the OPE 
is an operator equation by choosing the simplest possible states for the matching. 

The matrix elements of the operators in the OPE can be evaluated using diagrammatic methods by applying the 
Feynman rules given in Appendix A. The bilocal operator -0^.(1?— ^r)^p„{R+^r) can be represented diagrammatically 
by a pair of open dots with an atom line ending at the dot associated with tprr and an atom line beginning at the dot 
associated with ip^, as in Fig. 2(a). The atom lines should be labelled by the spin a, but we suppress that label. Local 
operators that annihilate one atom and create one atom, such as V-'J^cr and tpl'V^ipa — ^■''4't''Pcn can be represented 
diagrammatically by a single open dot, with an atom line ending at the dot and another atom line beginning at the 
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FIG. 3: Diagrammatic representations of the local operators (a) ^\^ltpii^2, (b) ^\4>l4i, (c) 4i^^l)i^p2, and (d) 




FIG. 4: The diagrams for matrix elements of the bilocal operator ijjl{R — ^r)tpa{R + ^r) between 2-atom scattering states. 
The open dots represent the operators. 

dot, as in Fig. 2(b). Local operators that annihilate a pair of atoms with spins 1 and 2 or a molecule and create a pair 

of atoms or a molecule, such as 'i/'|'i/'|V-'i'!/'2, ip\ip2't>7 4>^'>Pi'^-'2, and 4>^4>, can be represented diagrammatically by a single 
open dot, with two atom lines or a molecule line ending at the dot and two atom lines or a molecule line beginning 
at the dot, as in Fig. 3. 

The Wilson coefficients for operators that annihilate one atom and create one atom can be obtained by matching 
matrix elements between the 1-atom states consisting of a single atom with specified momentum and spin. For the 
bilocal operator 'tpl{R— ^r)ip„{R+ ^r), the matrix element between 1-atom states is given by the diagram in Fig. 2(a). 
The expression for the matrix clement is simply thc^ Feynman rule for the operator, which is an exponential function 
of r. For local operators that annihilate one atom and create one atom, such as ipl.ipa-(R), the matrix element is given 
by the diagram in Fig. 2(b). Their Wilson coefficients can be determined by matching the matrix elements between 
the 1-atom states. The Wilson coefficients are identical to those that would be obtained by multiplying the Taylor 
expansions of the operators ipl{R — ^r) and Vo-(-R + ^r). 

The Wilson coefficients of local operators that annihilate one atom and create one atom, such as V'J-V'cr (-R) ; could 
also be determined by matching matrix elements between 2-atom states, but the calculations are more complicated. 
The diagrams for the matrix element of the bilocal operator '(pl{R — ^r)^pcr{R + ^r) are shown in Fig. 4. The blobs 
represent the infinite series of diagrams that are summed up by the integral equation in Fig. 1. The contributions 
from each of the three diagrams in Fig. 4(a,b,c) are analytic functions of r. The diagram in Fig. 4(d) is evaluated in 
subsection A 2 a of the Appendix. It is an exponential function of r = |r|, so it is not analytic at r = 0. The diagrams 
for the matrix element of a local operator that annihilates one atom and creates one atom, such as tp'liJjrr{R), are 
shown in Fig. 5. Given the Wilson coefhcients determined by matching matrix elements between 1-atom states, the 
three diagrams in Fig. 5(a,b,c) match the expansions in powers of r of the corresponding diagrams in Fig. 4(a,b,c). 
The diagram in Fig. 5(d) matches the even powers of r in the expansion of the diagram in Fig. 4(d) in powers of 
r. The diagram in Fig. 4(d) is evaluated in subsection A 2 b of the Appendix. However for operators whose Wilson 
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FIG. 5: The diagrams for matrix elements of local operators that annihilate an atom and create an atom, such as V'ctV'ct or 
i>l^i>a- — ^i'ti'a, between 2-atom scattering states. The open dot represents the operator. 




FIG. 6: The lines in the diagram of Fig. 4(d) whose energies and momenta must be integrated over. The solid dots are vertices 
to which are attached either a pair of atom lines or a molecule line. 

coefRcient is an odd function of r, such as ip'^X/tpa — ^i^t'^Pa, the diagram in Fig. 4(d) vanishes. The odd powers of r 
in the expansion of the diagram in Fig. 4(d), which are not analytic at r = 0, must be matched by more comphcated 
operators whose matrix elements between 1-atom states arc zero. The next most complicated operators with nontrivial 
matrix elements between 2-atom scattering states have one factor of ipi'4'2 or and one factor of tpl'tl'l or 0^- 

To identify the operators on the right side of the OPE whose matrix elements between 1-atom states are zero but 
whose matrix elements between 2-atom scattering states arc nonzero, wc consider the lines in the diagram in Fig. 4(d) 
whose energies and momenta are integrated over. Those lines are shown in Fig. 6, where the momenta are labelled. 
The solid dots are vertices whose Feynman rule is —iXo/m if a pair of atom lines is attached and —igo/m if a molecule 
line is attached. Nonanalytic dependence on r can arise only from the integration region with large momentum q. 
There can be a significant contribution from this region only if the coordinates associated with the cndpoints of the 
lines carrying momentum q are close together. Contributions from this region therefore correspond to shrinking those 
lines to a single point. If we include the lines attached to the solid dots in Fig. 6, the shrinking of the large-momentum 
lines to a point gives one of the four vertices in Fig. 3. These four vertices correspond to the operators ipliplijjiijj2, 

(f)^iljiip2, V'lV'l'/'' ^^'^ '/'V) respectively. They appear with relative weights Ag, Aq^o, Aq^o, and g^, respectively. It is 
convenient to introduce a composite operator $(i?) defined by 

^R) = AoViV-f + goHR)- (43) 

The simplest combination of operators that corresponds to the large-momentum region of the diagram in Fig. 4(d) is 
therefore $1'$(il). 

The Wilson coefficient of the $t$(^) term can be determined by matching matrix elements between two-atom 
scattering states. The matrix element for $^$(ii) can be represented by the sum of the four diagrams in Fig. 7, 
together with 12 other diagrams in which there is no scattering of the two incoming lines or no scattering of the 
two outgoing lines or both. The sum of all these diagrams is calculated in subsection A 2 c of the Appendix. The 
dependence on the initial and final states is exactly what is required to match the term linear in r in the expansion 
of the diagram in Fig. 4(d). By matching, we determine the Wilson coefficient of the operator $1'$ to be — r/(87r). 
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FIG. 7: Diagrams for matrix elements of the local operator between 2-atom scattering states. The open dot represents 
the operator. For each of the diagrams shown, there are three other diagrams with no scattering of the atoms in the initial 
state or in the final state or both. 

Thus the analog of the OPE in Eq. (30) in the Resonance Model is 

ijtiR- ir)V.(i2+ ir) = i^lMR) + \r ■ {^l^^a{R) - VViV.(i2)] 

-^$t$(i2) + ..., (44) 

OTT 

Wc can now identify the coefficient C of the Xjk^ tail in the momentum distribution in Eq. (8). The leading term 
in the asymptotic expansion of the momentum distribution in Eq. (29) comes from the term with operator $t<I>. The 
Fourier integral of its Wilson coefficient — r/(87r) is given by Eq. (31). Comparing with Eq. (8), we find that the 
contact in the Resonance Model is 

G = Jd^R{X\^^^R)\X), (45) 

where $(i2) is the composite operator defined by Eq. (43). In subsequent sections, we will suppress the state \X) and 
denote the expectation value of an operator in that state simply by {O). 



VI. ENERGY RELATION 



Tan's energy relation in Eq. (9) is a decomposition of the energy into terms that are separately insensitive to the 
range of the interactions between the atoms. In a quantum field theoretic formulation of the problem, the analogous 
relation is a decomposition of the energy into terms that are insensitive to the ultraviolet cutoff. To accomplish 
this separation for the Resonance Model, we use the facts that and are operators with ultraviolet-finite 
matrix elements. For this follows from it being the contact density operator. For this follows from an 

adiabatic relation derived later in Section VIII. Using the expressions for the bare coupling constants in Eqs. (38), 
the interaction term T in the Hamiltonian density in Eq. (37b) can be expressed in the form 

Z= fi/- - -4—^ (46) 

We can now decompose the sum of the kinetic energy density and the interaction energy density into five terms, each 
of which has ultraviolet-finite matrix elements: 

^ + ^ = i E (v^i • - ^*^*) + i v^t . 



(47) 
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To sec that the term that is summed over a in Eq. (47) has ultraviolet finite matrix elements, we consider the integral 
over space of its expectation value: 



(48) 



This differs from the expectation value of the kinetic energy of atoms with spin cr by a subtraction that is hncar in the 
ultraviolet cutoff A. By expressing the quantum fields ij: and i/)''' in terms of their Fourier transforms and using the 
expression for the momentum distribution Pcr(fc) in Eq. (26), we can express the subtracted kinetic energy in Eq. (48) 
in the form 



(49) 



where C is the contact given in Eq. (45). We have expressed the factor of A in Eq. (48) as an integral over momentum 
space. Since the tail of the momentum distributions at large k has the form in Eq. (8), the subtraction in Eq. (49) 
makes the integral over k convergent in the limit A ^ oo. 

The energy relation for the Resonance Model is obtained by taking the expectation value ofW = T + J+ V and 
integrating over space. It can be expressed in the form 



where {V) is the energy associated with the external potential, 



mX 



{V)= d'RViR){4i;i+4i;2 + 2cP^), 



and (Tmoi) is the kinetic energy of the molecules: 



(50) 



(51) 



(52) 



VII. LOCAL PAIR DENSITY 



The Tan relation in Eq. (10) expresses the contact density C in terms of the local pair density Npa,ij.{R, s), which is 
the number of pairs of atoms in a ball of volume |7rs^ centered at R. To obtain the analogous relation in the Resonance 

Model, we consider the OPE of the number density operators 4tjji{R — |r) and ip2'^2iR+ ^r). The matrix elements 
of the product of these operators between 2-atom scattering states can be represented diagramatically by the sum of 
the four diagrams in Fig. 8. The three diagrams in Figs. 8(b,c,d) involve integrals over the 4-momenta of atoms. The 
integrated momenta run through the lines that connect the operator vertices to the blobs. Wilson coefficients that are 
not analytic at r = must come from the large-momentum regions of those integrals. The corresponding operators 
can be deduced by shrinking the lines carrying the large momentum to a point. Shrinking those lines in the diagram 
in Fig. 8(d) produces the operator vertices in Fig. 3 in the combination that corresponds to the operator $^<S(il). 
The Wilson coefficient of this operator is calculated in subsection A2d of the Appendix. It is proportional to l/r^, so 
it is singular at r = 0. For the diagrams in Fig. 8(b) and (c), shrinking the lines carrying large momentum produces 
operator vertices that correspond to the operators ipltp2^{R) and $'l'^i?/;2(-R), respectively. The Wilson coefficient of 
these operators are proportional to 1/r. Thus the most singular term in the OPE is 

We can use the OPE in Eq. (53) to obtain a simple interpretation of the contact density C. It can be expressed in 
the form 

where on the right side we have kept only the most singular term as r' r. If we integrate the left side of Eq. (54) 
over both r and r' inside the ball of volume Itts"^ centered at R, we obtain an operator that in the absence of 
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(a) 




(c) 




id) 




FIG. 8: The diagrams for the matrix clement of the bilocal operator ■tp\tpi{R — ^r)ipl'4>2{R + ^r) between 2-atom scattering 
states. The open dots represent the operators. 

interactions would give the product of the number of atoms with spin 1 and the number of atoms with spin 2 in that 
ball. Thus the expectation value Afpair(^i, s) of this operator can be interpreted as the number of pairs of atoms with 
spins 1 and 2 inside that ball. Integrating the right side of Eq. (54) over r and r' inside that same ball and taking 
the expectation value, we obtain the limiting behavior of A^pair(^, s) as s — > 0: 



7Vpair(-R,s) 



(55) 



Since (<i>^$(il)) is the contact density C, the Tan relation in Eq. (10) is satisfied without modification in the Resonance 
Model. 



VIII. ADIABATIC RELATIONS 



Tan's adiabatic relation in Eq. (11) gives the rate of change in the energy £^ of a state due to a change in the 

interaction parameter a. In the Resonance Model, there are three interaction parameters: ly, g, and A. According to 
the expressions for these parameters in Eq. (36), v can be changed experimentally by varying the magnetic field while 
g and A are essentially constants. Nevertheless, we will consider the effects of changing all three parameters. 
The Feynman-Hellman theorem for variations in the parameter v is 



dE/du = jd^R{dn/diy). 



(56) 



There are analogous equations for the variations in the parameters g and A. Since the Resonance Model is renormal- 
izable, the Hamiltonian density operator H = T + T + V given by Eqs. (37) has ultraviolet-finite matrix elements. 
More explicitly, if Eqs. (38) are used to eliminate the bare parameters vq, g^, and Aq in favor of the renormalized 
parameters v, g, and A, the matrix elements of 7i have finite limits as A — » oo, where A is the ultraviolet cutoff. Since 
H has ultraviolet-finite matrix elements for any values of u, g, and A, its derivatives {d/du)'H, {d/dg)H, and {d/dX)'H 
must also have ultraviolet-finite matrix elements. The partial derivatives with the ultraviolet cutoff A held fixed can 
be evaluated by using Eqs. (38) for the bare parameters. The expressions for the partial derivatives are simplest if 
they are expressed in terms of the field (p and the composite operator $ defined in Eq. (43): 



d 



2.9^ 



(57a) 
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The right sides have been expressed as hnear combinations of the operators (t)^<i), + and ^^cj)^ with coefficients 
that are functions of renormahzed parameters. Since the three partial derivatives of TL in Eqs. (57) must have 
ultraviolet-finite matrix elements, these three operators must also have ultraviolet-finite matrix elements. 

The adiabatic relations for the Resonance Model can be obtained by taking the expectation values of the partial 
derivatives in Eqs. (57) and then integrating over all space: 

y^E = V J (f R (j)), (58a) 



IX. PRESSURE RELATION 



For a homogeneous system, the external potential is a constant: V{r) = — /x, where fi is the chemical potential. 

The thermodynamic properties of the system are determined by the free energy density (or thermodynamic potential 
density) It can be expressed in terms of the trace of the partition function: 



exp(-/?JPy) = Trexpi-P J d^RH), 



(59) 



where (5 =1/T and V is the volume. The free energy density is a function of T, /x, and the interaction parameters 
V, g, and A. The parameter u has dimensions of energy. The mass m, together with h = 1, can be used as a conversion 
constant to obtain energy scales l/(mA^) and g'^/in associated with the two other interaction parameters. The free 
energy density has dimensions of energy per volume, which can be expressed as (energy) ^/-^ up to factors of m and 
U. Thus dimensional analysis requires that the free energy density satisfy 



d 



d 



d 



Id IB 



T = -T 
2 



(60) 



The logarithmic derivatives on the left side of Eq. (60) can be determined by applying them to the logarithms of 
both sides of Eq. (59). In the case of the interaction parameters i^, g, and A, the logarithmic derivatives of T are 
simply the expectation values of the logarithmic derivatives of the Hamiltonian density given in Eqs. (57). In the case 
of the temperature and the chemical potential, the logarithmic derivatives of T are 



dT 
d 



— £ + jin, 



where n is the expectation value of the total number density: 

n = (V'lV'i + -02^2 + 2 



(61a) 
(61b) 

(62) 



Because the Hamiltonian density Ti includes the term V in Eq. (37c), it is actually the grand Hamiltonian density 
for the grand canonical ensemble. Its expectation value is {H) = £ — /in, where £ is the energy density. Because the 
system is homogeneous, the pressure is simply V = —T. Combining together all the terms in Eq. (60), we obtain 
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1 



4mA ^' 2mX 



(63) 



This is the pressure relation for the Resonance Model. 



X. VIRIAL THEOREM 



If the trapping potential V{R) has simple behavior under scaling the coordinates R = {X, Y, Z), the virial theorem 
gives simple relations between various contributions to the energy. The most important case for cold atoms is a 
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harmonic trapping potential. For simplicity, we consider only the case of an isotropic harmonic potential: 



1 



V{R) = -muj^R\ 



(64) 



The generalization to the case of an anisotropic harmonic potential is straightforward. 

The virial theorem can be derived simply from dimensional analysis. Let us assume that the state of the system can 

be specified by a set of integer quantum numbers. The energy E = J (Pr{T-l) for any such state must be a fimction of 
the frequency uj and the interaction parameters u, g, and A. Since E has dimensions of energy, dimensional analysis 
requires that it satisfy 
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E = E. 



(65) 



If we consider the system at a nonzero temperature T, the differential operator in Eq. (65) would also include a term 

T{d/dT). 

The action of the differential operator in Eq. (65) on the Hamiltonian density can be obtained by inserting the 
partial derivatives in Eqs. (57): 



' d d Id 15 



9 



n = 2V{R) (v-iVi + ^2V'2 + 2(j)U 



mX 



3g 

4mA 



(<^t$ + $t0) 



2toA 



(66) 



The virial theorem for the Resonance Model is obtained by taking the expectation value, integrating over space, and 
setting this equal to E: 



E 



= 2{V)+ jSR 



9 



35 



(67) 



where {V) is the energy from the external potential given by Eq. (51). 



XI. INELASTIC TWO-BODY LOSS RATE 



If the atoms have other spin states with lower energy than the spin states a = 1, 2 of interest, a pair of atoms with 
spins 1 and 2 may be able to scatter inelastically into the lower spin states. The molecular state responsible for the 
Feshbach resonance may also be able to decay into pairs of atoms in lower spin states. If the lower spin states interact 
weakly with the spin states 1 and 2, the primary effects of transitions to the lower spin states is the disappearance of 
atoms from the spin states 1 and 2 and the disappearance of the Feshbach molecule. The inclusive effects of the lower 
spin states on states near the threshold for the spin states 1 and 2 can be taken into account through anti-hermitian 
terms in the Hamiltonian density. These terms must be anti-hermitian to account for the loss of probability from 
transitions to the lower spin states, which are not taken into account as explicit degrees of freedom. 

If the energy gap between the lower spin states and the spin states 1 and 2 is much larger than the energies 
of interest, the anti-hermitian terms in the Hamiltonian density can be chosen to be local operators. The rate of 
disappearance of molecules through their decays into pairs of atoms in lower spin states is proportional to the number 
density of the molecules. The resulting effects can be taken into account through the operator (j)^(f> with a negative 
imaginary coefficient. The rate of disappearance of pairs of atoms through inelastic scattering is proportional to the 
local density of pairs of atoms. The resulting effects can be taken into accoimt through the operator with a 

negative imaginary coefficient. If the decay of the molecule and the inelastic scattering of a pair of atoms can produce 
the same final states, the resulting interference effects can be taken into account through the operators (f)^ip\4'2 and 
tpltplcf). The large energy gap -Egap is essential for allowing the effects of transitions into the lower spin states to 
be taken into account through local operators. The atoms created by the decay of a molecule or by the inelastic 
scattering of atoms with spins 1 and 2 will emerge with large momentum given approximately by (mSgap)^/^. By the 
uncertainty principle, a localized wave packet corresponding to such an atom can be traced back to a region of origin 
whose size is approximately Ti/ {mE^^^Y^^ . But the regions of origin of the two atoms must also be separated by a 
distance at most of order h,/ {mEg-^^y/^. The reason for this is that the interactions that change the spin state of the 
atom from 1 or 2 to a lower spin state must also deliver a momentum kick of order [mEga^pY^'^ . This large momentum 
can be transferred between the two final state atoms only if they are close together. Even if it is transferred though 
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an intermediate atom with spin 1 or 2, that atom must be far off its energy shell and therefore can propagate only 
over a short distance. As long as %/ {mE^^pY^'^ is small compared to the length scales that are described explicitly 
by the field theory, the effects of the transitions to lower spin states can be taken into account systematically through 
local operators. 

The operators V'iV'2V'iV'27 + V'lV'l'^? (f)^ (j) required to take into account the leading effects of the lower 

spin states are ones that appear already in the interaction term in the Hamiltonian density given in Eq. (37b). To 
give anti-hermitian terms in the Hamiltonian density, these operators must be multiplied by imaginary coefficients. 
Thus the leading effects of the lower spin states can be taken into account by adding negative imaginary parts to 
the bare parameters z/o , .go , and Aq . We now consider the effect of the nonhermitian terms in the Hamiltonian. The 
expectation value of those terms can be expressed as — ir/2, where 



-2 jd^R [lm..o(0V) + ^(0Vi^f)+V'IV'|0(^^) + ^(^I4V'i^f^) 



(68) 



The leading effect of these terms on a state with definite energy E is to change its time-dependence from cxp{—iEt/h) 
to exp(— z(i? — iT/2)t/h). The probability in such a state decreases with time at the rate T/h. Thus F is the rate at 
which atoms with spins 1 and 2 and Feshbach molecules are depleted by transitions to pairs of atoms in lower spin 
states. 

The loss rate in Eq. (68) can be expressed in terms of expectation values of operators with ultraviolet-finite matrix 
elements by eliminating 'ipiip2 in favor of the operator defined in Eq. (43). Since the coefficients of these expectation 
values must also be ultraviolet finite, they are expressible in terms of the renormalized parameters v, g, and A. Using 
Eqs. (38) to eliminate the bare parameters and expanding to first order in the imaginary parts of the renormalized 
parameters, we obtain 



+ ( ^Im.g ^ImA ) <^ 



V niX 



iA2 



ImA($T$\ 



(69) 



This is the generalization in the Resonance Model of the Tan relation for the inelastic 2-body loss rate. 



XII. CONSTRAINTS FOR LARGE SCATTERING LENGTH 



The Tan relations in Eqs. (8)-(14) express various properties of an arbitrary state consisting of atoms with two spin 
states and a large scattering length in terms of the contact C or the contact density C. These relations hold exactly 
in the Zero-Range Model. The contact density operator in this model was identified in Ref. [10] to be AqViV'2V'iV'2- 
In the Resonance Model, the generalizations of the Tan relations derived in Sections V-XI involve the expectation 
values of three distinct local operators: $t<j)_ (^1$ _|_ $t0^ and i/it^, where $ is the composite operator defined in 
Eq. (43). The contact density operator is identified in Eq. (45) to be The same operator appears in the energy 

relation in Eq. (50) and in the local pair density in Eq. (55). The three adiabatic relations in Eqs. (58), the virial 
theorem in Eq. (67), the pressure relation in Eq. (63), and the 2-body inelastic loss rate in Eq. (69) involve various 
linear combinations of the expectation values of the three local operators. 

The scattering length is the only interaction parameter in the Zero-Range Model. The three parameters of the 
Resonance Model provide three distinct length scales: \mi'\~^/^, g~'^ and |A|. There are two limits in which the 
scattering length given by Eq. (35a) becomes arbitrarily large compared to the other two length scales. In particular, 
it becomes large compared to the effective range given by Eq. (35b). One of these limits is A ^ ±oo with v and g 
fixed, in which case a — > A/(47r). The other limit is — > with g and A fixed, in which case a —g'^/{A'Kmu). In 
either of these two limits, the universality of systems with large scattering length requires the Resonance^ Model to 
reduce to the Zero-Range Model. The generalized Tan relations must therefore reduce to the original Tan relations. 
We will find that this imposes constraints on matrix elements of and -|- $^0. 



A. A ^ ±oo with 1/ and g fixed 



We first consider the case of a large scattering length a obtained by increasing | A| with v and g fixed. The adiabatic 
relation in Eq. (58c) can be expressed as an adiabatic relation for variations in a: 

j d^R {g^<j>U) - 9{<l>^^ + + (^^^)) • (70) 
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When A is sufficiently large, this must reduce to Tan's adiabatic relation in Eq. (11), which can be expressed as 

Le=^ [ d^R ($t$) + 0(1/a3). (71) 



da mX^ 

This implies a constraint on the matrix elements in Eq. (70) as A — > ±c»: 
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2uK 



. g{(l)-^^ + <^'<(l)) = 0{1/X). (72) 

We can obtain an independent constraint on the matrix elements by using the virial theorem, which for the Resonance 
Model is given in Eq. (67). When A is sufficiently large, this must reduce to Tan's virial theorem in Eq. (13), which 
can be expressed as 

E = 2{V)--^ /d3i?($t$)+0(i/A2). (73) 
This implies a constraint on the matrix elements in Eq. (67) as A ^ ±oo: 

^ " ^'^^'^^ + 4^^^^*^ + '^^'^^ = ^^^^ 

If 1/ 7^ 0, the constraints in Eqs. (72) and (74) imply 

(0t0) = 0(1/A2), (75a) 
(0t$ + $t^^ = 0(1/A). (75b) 

These constraints are consistent with the decoupling of the molecular field (f) in the limit A — > ±oo. Matrix elements 
of operators involving the field (j) are suppressed by a factor of 1/A for every factor of (f) in the operator. 



B. 1/ — > with g and A fixed 

We next consider the case of a large scattering length a obtained by decreasing \v\ with g and A fixed. According 

to Eq. (36a), this can be achieved by tuning the magnetic field to the position of the Feshbach resonance. The 
adiabatic relation in Eqs. (58a) can be expressed as an adiabatic relation for variations in a: 



d Airmv^ 
-E 



J d^R {(!>'' (/)). (76) 



da 

When V is sufficiently small, this must reduce to Tan's adiabatic relation in Eq. (11), which can be expressed as 



d Airmu'^ 
-hi = 



jd^R{<^^^) + 0{i^^). (77) 



da 

This implies that in the limit — > 0, the expectation value of <f)^<f) can be expressed in terms of the contact density: 

g2(0V) = (ci>t$)+o(^). (78) 

We can obtain an independent constraint on the matrix elements by using the virial theorem, which for the 
Resonance Model is given in Eq. (67). When u is sufficiently small, this must reduce to Tan's virial theorem in 
Eq. (13), which can be expressed as 

E = 2{V) + ^ Jd^R + 0(i.2). (79) 

This implies a constraint on the matrix elements in Eq. (67) as ^ 0: 

Combined with the constraint in Eq. (78), this implies that the expectation value of + can also be expressed 
in terms of the contact density: 

^(^■f^ + ^V) =2($^$) + 0(j^). (81) 
It also implies a constraint on the terms of order i/ on the right sides of Eqs. (78) and (81): 

{g\cl>U) - (*^*)) - I (5(0^^ + ^U) - 2($^$)) = ^(^^^) + 0{u'). (82) 
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XIII. EFFECTIVE RANGE MODEL 



If we set A = in the S-wave phase shift for the Resonance Model, which is given by Eq. (6), it reduces to 

5^ 9 



fccot(5o(A;) = — ^2 2^^- (83) 



This is the phase shift of the Effective Range Model in Eq. (7), with the scattering length and the effective range 
given by 

r, = (84b) 

Note that the effective range is negative. 

In the simplest quantum field theory formulation of the Effective Range Model [20], there are three quantum fields: 
fermionic fields ipi and il'2 that annihilate atoms and a bosonic field (p that annihilates a diatomic molecule. The 
Hamiltonian density 7i = T + 1 + V for the Effective Range Model is the sum of the kinetic term in Eq. (37a), the 
external potential term in Eq. (37c), and the interaction term 

I = uo<l>U + ^ ('/'ViV'f^ + i^UU^^^) ■ (85) 

To avoid ultraviolet divergences, an ultraviolet cutoff A must be imposed on the momenta of virtual particles. The 
only renormalization required in this model is an additive renormalization of the parameter uq: 

uo{A) = ,.+ ^. (86) 

The Effective Range Model can also be formulated as a local quantum field theory with fermion fields tpi and tl>2 only 
and with a momentum-dependent interaction [20]. In this formulation of the model, a rather intricate renormalization 
of the parameters is required to obtain the phase shift in Eq. (83) in the limit A ^ oo. 

The generalized Tan relations for the Effective Range Model can be derived by using the same methods that we used 
for the Resonance Model. However they can be derived more easily by starting from the generalized Tan relations 
for the Resonance Model and taking the limit A ^ 0. This limit is not trivial, because some of the generalized Tan 
relations for the Resonance Model involve factors of 1/ A. The limit A ^ can be taken by expanding the composite 
operator $ defined in Eq. (43) through first order in A. Using the expressions for the bare coupling constants in 
Eqs. (38a) and (38b) to expand Aq and go in powers of A, we obtain 



$(R) = #(R) + A 

We denote the composite operator defined by the coefficient of A by $'(i?): 



0(A2). (87) 



^'{R) = ViVf ^(i2) + ^m)- (88) 

The contact C in the Resonance Model is given in Eq. (45). From the expansion in Eq. (87), we can see that the 
limit A — > can be taken simply by replacing $ by g(p. Thus the contact in the Effective Range Model is 

C = j SRg'^{X\(l>'^(l>{R)\X). (89) 

We proceed to enumerate the generalized Tan relations for the Effective Range Model: 

1. Tail of the momentum distribution. This is given by Eq. (8), where C is the contact in Eq. (89). 

2. Energy relation. The decomposition of the energy into terms that are ultraviolet-finite can be obtained by 
taking the limit A ^ in Eq. (50): 

E = {V) + (Ti^^"^') + (Ti^"^)) + (T^oi) + j d^R {y{<P^<t>) + ^ (<^^*' + V)) , (90) 
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where <J>' is the composite operator defined by Eq. (88). The subtracted kinetic energy for a single spin state is 
obtained from Eq. (49) by replacing $ by 

^y(sub)^ = ^Jd'R (^(Wj • VV^S^)) - • (91) 

3. Local pair density. The contact density C = g'^{(j)^4') is related to the local pair density A^pair(-R, s) by Eq. (10). 

4. Adiabatic relation. The adiabatic relation for variations in v is the same as in the Resonance Model, where 
it is given by Eq. (58a). The adiabatic relation for variations in g can be obtained from Eq. (58b) by taking the 
limit A — > 0. The adiabatic relations are therefore 

u^E = V jd^R{4>''(l)), (92a) 

g^E = ^ fd^R{<f>^^' + ^''i<i)). (92b) 
dg m J 



5. Pressure relation. The relation between the energy density £ and the pressure "P in a homogeneous system 
can be obtained from Eq. (63) by taking the limit A — > 0: 

f = ^p + :.(0t0) + _|.(0t$' + $'t^). (93) 

6. Virial theorem. The virial theorem for a system with a harmonic trapping potential can be obtained from 
Eq. (67) by taking the limit A 0: 

E = 2{V)+ jd^R(i^{(l)U) + £^{<t>^^' + ^'U))- (94) 

The virial theorem for the Effective Range Model has been derived previously by Werner [26]. His result can be 
expressed in the form 

This is just the generalization of the dimensional analysis condition for the energy in Eq. (25) for a model with 

two interaction parameters a and with dimensions of length. If we use Eqs. (84) to change variables from a 
and to v and g and if we use the adiabatic relations in Eqs. (92) to express the derivatives in terms of matrix 
elements, we obtain the virial theorem in Eq. (94). 

7. Inelastic 2-body loss rate. The rate at which low-energy atoms and molecules are depleted by transitions to 

pairs of atoms in lower spin states can be obtained from Eq. (69) by setting ImA = and then taking the limit 
A ^ 0. It can be obtained more simply by using the adiabatic relations in Eqs. (92) to determine the change in 
the energy due to small imaginary changes in the parameters u and g: 



r«-2 Jd^R (^lmiy(^V) + ^^(<^^$' + $'^ 



(96) 



The scattering length in Eq. (84a) can be made arbitrarily large by taking the limit with g fixed. In this 
limit, the universality of systems with large scattering length requires the Effective Range Model to reduce to the 
Zero-Range Model. This imposes a constraint on the expectation values of the two operators that appear in the 
generalized Tan relations. The constraint can be derived by demanding that the virial theorem in Eq. (94) reduce to 
the Tan relation in Eq. (13) as ^ 0. It can also be derived from the constraint for the Resonance Model in Eq. (80) 
by taking the limit A ^ 0. They both give the same constraint: 



(97) 
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XIV. SUMMARY 

The Tan relations enumerated in Section II are universal relations between various properties of an arbitrary system 

consisting of fcrmions with two spin states and a large scattering length. These relations should be satisfied provided 
the energy, temperature, and number density of the system are low enough that the interactions between the atoms 
are accurately described by the S-wave phase shift in Eq. (5), which depends only on the scattering length a. They are 
satisfied exactly in the Zero-Range Model, in which the phase shift is given by Eq. (5) up to arbitrarily high energies. 

In Ref. [10], quantum field theory methods were used to derive the Tan relations. The Zero- Range Model can be 
formulated as a local quantum field theory. Using this formulation, the Tan relations were derived by using standard 
renormalization methods together with the OPE. The contact density operator in this model was identified in Ref. [10] 
to be AqV'i^2V'iV'2- One advantage of using quantum field theory methods is that it is straightforward to derive the 
generalizations of the Tan relations for any system that can be formulated as a renormalizable local quantum field 
theory. In this paper, we used quantum field theory methods to derive the generalized Tan relations for the Resonance 
Model. We also derived them for the Effective Range Model by using the fact that it can be obtained as a limit of 
the Resonance Model. 

The Resonance Model is defined by the S-wave phase shift in Eq. (6), which depends on three interaction parameters 
u, g, and A. In Sections V-XI, we derived the generalized Tan relations for the Resonance Model. We identified the 
contact density operator to be $t<I>, where $ is the composite operator defined in Eq. (43). The generalized Tan 
relations involve expectation values of linear combinations of three local operators: cp^^ + 4>^4>- The 

scattering length in the Resonance Model can be made arbitrarily large by tuning i/ — > 0. The condition that the 
generalized Tan relations reduce in this limit to the original Tan relations implies the constraints on the expectation 
values of <i>^<i>, (j)^^ + ^^cj), and that are given in Eqs. (78), (81), and (82). 

The Effective Range Model is defined by the S-wave phase shift in Eq. (7), which depends on the two parameters 
a and rg < 0. It can be obtained from the Resonance Model by setting A = 0, in which case the S-wave phase 
shift is given by Eq. (83) in terms of the parameters u and g. In Section XIII, we enumerated the generalized Tan 
relations for the Effective Range Model. We identified the contact density operator to be g'^(f>'^(f>. The generalized 
Tan relations involve expectation values of linear combinations of two local operators: (p'^cp and (p^^' + ^ ^(p, where 
is the composite operator defined in Eq. (88). The scattering length in the Effective Range Model can be made 
arbitrarily large by timing — ^ 0, which implies \a\ ^ [r,,!. The condition that the generalized Tan relations reduce 
in this limit to the original Tan relations implies the constraint on the expectation values of ^'^^ and (j)^^' + ^ ^(j) that 
is given in Eq. (97). 

The Tan relations can be tested experimentally by using cold trapped atoms near a Feshbach resonance, as discussed 
at the end of Section II. Measurements of the momentum distributions can be used to determine the contact C and the 
sum (T) + (I) of the kinetic and interaction energies. Measurements of the density profiles can be used to determine 
the potential energy {V). The virial theorem in Eq. (13) is a nontrivial relation between (T) + (/), {V), and C. 
The adiabatic relation in Eq. (11) gives a constraint on the variations in (T) + (/), (V^), and C from changing the 
scattering length. Both of these relations can be tested experimentally. Since the ability to change the scattering 
length is essential both for determining C and (T) + (/) and for testing the adiabatic relation, the Feshbach resonance 
plays a crucial role. The Resonance Model provides a more detailed microscopic model for atoms near a Feschbach 
resonance than the; Zero-Range Model. Our results on the generalized Tan relations in the Resonance Model should 
be useful for quantifying the theoretical errors in experimental tests of the Tan relations. 
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APPENDIX A: DIAGRAMMATIC CALCULATIONS 
1. Feynman rules 

The results for the matrix elements given in Sections V and VII can be calculated relatively easily by applying a set 

of rules - known generally in quantum fic^ld theory as Feynman rules. The calculation of matrix elements of operators 
reduces to drawing all relevant Feynman diagrams, using the Feynman rules to write down mathematical expressions 
for the matrix elements, and evaluating any integrals that appear in those expressions. 
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a. Propagators and vertices 

We first give the Feynman rules for tlie Resonance Model that are required to calculate Green's functions and 

T-matrix elements. A Feynman diagram consists of atom lines (represented by single lines) and molecule lines 
(represented by double lines) connected by vertices of two types. The mathematical expression for a diagram is 
obtained by applying the following rules: 

1. Four-momenta. Assign a four- momentum {po,p), where is the energy and p is the momentum, to every 

external line. The initial and final four-momenta are constrained by overall conservation of energy and momen- 
tum. The four-momentum of the internal lines are constrained only by conservation of energy and momentum 
at each vertex. 

2. Propagator factors. For each internal line with four-momentum [q^, q), include the factor i/[qQ — (^ /{2m) + ie\ 
for an atom line and i/[qo — v — /{Am) + ie\ for a molecule line. 

3. Vertex factors. For each 2-atom-to-2-atom vertex, include the factor —iXo/m. For each 2-atom-to-molecule 
vertex, include the factor —igo/m. 

4. Loop momenta. If there are four- momenta {ko,k) that are not determined by the four-momenta of the 
external lines, integrate over them using the measure / kdko / {2ttY . The integrals over ko can be evaluated 
using contour integration. 

The connected Green's function for N atoms to evolve into N atoms can be expressed as the sum of all connected 
diagrams with N incoming atom lines and outgoing atom lines. An example is the amplitude A{E) for AT = 2 
in Eq. (42), which can be expressed as the infinite sum of diagrams generated by iterating the integral equation in 
Fig. 1. If this amplitude occurs as a subdiagram of a Feynman diagram, the Feynman rule for the subdiagram is 
iA{E — P^/(4m)), where E and P are the total energy and momentum flowing through the subdiagram. 

The T-matrix element for N atoms to scatter into N atoms is obtained by setting the energy for each external 
atom line equal to the value required by the nonrelativistic energy-momentum relation: po = p'^ /(2m). The T-matrix 
elements involving diatomic molecules are more complicated and will not be discussed here. 



b. Operator vertices 



To calculate matrix elements of composite operators, we also need Feynman rules for the operators. Energy and 
momentum can flow into and out of the vertices for the operators. For simplicity, we consider only operators at a 
fixed time t = 0. The operators listed below in Eqs. (Al) annihilate an atom and create an atom, so the operator 
vertex has an incoming atom line and an outgoing atom line. If the incoming and outgoing momenta are k and k', 
the Feynman rules for the operator vertices are 



i;l{R-lr)MR+^r) 



exp[zfc • [R + ir)] exp[-ife' • {R - \r)], 

exp[i(fc - k') ■ R], 

i{k + ky exp[i{k - k') ■ R]. 



(Ala) 
(Alb) 
(Ale) 



The operators listed below in Eqs. (A2) annihilate a pair of atoms or a molecule and create a pair of atoms or a 
molecule. If the total momenta entering and leaving the vertex are K and K', respectively, the Feynman rules for 
the operator vertices are 



44iP,MR) 
(P^i^iMR) 

0t0(i2) 



exp{i{K — 


K') 


■R), 


exp{i{K — 


K') 


■R), 


exp{i{K — 


K') 


■R), 


exp{i{K — 


K') 


R). 



(A2a) 

(A2b) 
(A2c) 
(A2d) 



We will consider matrix elements of these operators only between states for which the total momenta K and K' are 
both zero, in which case the Feynman rules in Eqs. (A2) all reduce to 1. 

The matrix element of a composite operator between initial and final states consisting of atoms is the sum of all 
diagrams in which each external atom line is connected either to an operator vertex or to other external atom lines. 
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2. Matrix elements between 2-atom scattering states 

In this subsection, we calculate the matrix elements of composite operators between 2-atom scattering states. We 
denote the state consisting of a single atom with momentum p and spin a by |p, cr). A 2-atom scattering state is 
labelled by the momenta and the spins of the two particles: \pi,o'i;P2,'^2)- For simplicity, we will calculate matrix 
elements only between 2-atom scattering states consisting of two atoms with different spins, total momentum 0, 
and total energy E. These states, which we denote by |p, 1;— p, 2), are labelled by a vector p whose magnitude is 
\p\ = {mEy/^. 

a. Matrix element of tl'l{R — ^r)il>tT{R+ ir) 

The matrix element of the bifocal operator ip^{R — ^r)tpa{R + ^r) is given by the sum of the four diagrams in 
Fig. 4. The only one of these diagrams that involves an integral over the 4-momentum of an atom is Fig. 4(d). Using 
the Feynman rules, the expression for this diagram is 

ip', 1; -p', 2| - lr)MR + lr)\p, 1; -p, 2) ^^^^ 

= j2/pN / Arf^o expjiq ■ r) 

^ V (2^)4 [qo - qy{2m) + le^E - qo - gV(2m) + ie]' ^ ^ 

where A{E) is the amplitude in Eq. (42) and E — p'^/m. The integral over qo can be evaluated using contours. The 
resulting momentum integral is given in Eq. (A13b). The final result is 

2 

ip', 1; -p', 2\^UR - \r)MR + \r)\p, 1; -p, 2) |^^^^ = "^A^E) exp(»pr)- (A4) 
Unlike the first three diagrams in Fig. 4, this diagram is not an analytic function of r at r = 0. 

h. Matrix element of iplip^ 

The matrix element of the local operator il)^'ip^{R) is given by the sum of the four diagrams in Fig. 5. These 

diagrams are equal to the r — > limits of the corresponding diagrams for the bilocal operator in Fig. 4. For the 
first three diagrams, this can be seen easily from the Feynman rules. This is not as obvious for the fourth diagram, 
because it involves an integral over the 4-momentum of an atom. Using the Feynman rules, the expression for the 
diagram in Fig. 4(d) is 

ip', 1; -p', 2\i;iMR)\P, 1; -P, 2)|5(rf) 

>J (27r)4 [?o-gV(2m) + *e]2[i?-go-gV(2m) + ze]- ^ 

The integral over q^ can be evaluated using contours. The resulting momentum integral is given in Eq. (A14b). The 
final result is 

(pM;-p',2|VI^.(H)|p,1;-p,2)|^(,) = '-^A\E). (A6) 

This result matches the r° term in the expansion of Eq. (A4) in powers of r. This is consistent with the Wilson 
coefficient of il)^ip„{R) in the OPE in Eq. (44) being simply 1. 

c. Matrix element of ^ 

The matrix element for $^$(fi) can be represented by the sum of the four diagrams in Fig. 7, together with 12 
other diagrams in which there is no scattering of the two incoming atoms or the two outgoing atoms or both. Using 
the Feynman rules, the expression for the diagram in Fig. 7(a) is 

ip', 1; -p', 2| AgvIV'IV'iV'2(i?)|p, 1; -P, 2) I 



■^0 



■iAiE) 



7(a) 

(fiqdqo 1 



(27r)4 [go - g7(2m) + ie][E - qo - q^/{2m) + ie] 



(A7) 
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The integral over qo can be evaluated using contours. The resulting momentum integral is given in Eq. (A14a). The 

result is 



(p',l;-p',2|A2^IV|ViV'2(H)|p,l;-p,2)|7( . = 



mA{E) 



A ip 



(A8) 



The expressions for the three diagrams obtained from Fig. 7(a) by omitting the scattering of the atoms in the initial 
state or final state or both arc obtained from Eq. (A8) by replacing one or both of the factors in square brackets by 
1. The expressions for the three diagrams in Figs. 7(b-d) are obtained by replacing one or both of the factors of Aq 
by {gQ/m)/{E — uq). The sum of all 16 diagrams is 



{p', 1; -p', 2|$t$(ij)|p, 1; 2) = ( Ao + 



2 r 



1 + mA{E) 



By using the integral equation in Eq. (41), this can be simplified to 

(p', 1; -p', 2|$t<i,(i2)|p, 1; 2) = m'A^{E), 



A ip 



(A9) 



(AlO) 



where A{E) is the amplitude in Eq. (42) and E = p^ jm. This result matches the term in the expansion of Eq. (A4) 
in powers of r if the Wilson coefficient of $t$(i?) in the OPE in Eq. (30) is -r/(87r). 



d. Matrix element of tplipi{R — ^r)'4>2ip2{R + |r) 



The matrix element of the bilocal operator iplipi{R— ^r)ipltfj2{R + ^r) is given by the sum of the four diagrams in 
Fig. 8. The diagram in Fig. 8(a) is simply 1. The two diagrams in Fig. 8(b,c) involve an integral over the 4-momentum 
of an atom.. The diagram in Fig. 8(d) involves two such integrals. Using the Feynman rules, the expression for this 
diagram is 



(p', 1; -p', 2|vIVi(i? - kr)i'lMR + 1; p, 2) 

cfqdqo exp(iq ■ r) 



8(d) 



-iA{E) 



(2^)4 [qo - g2/(2m) + le] [E - qo - g2/(2m) + le] 



(All) 



The integral over qo can be evaluated using contours. The resulting momentum integral is given in Eq. (Al3a). The 
final result is 



(p', 1; -p', 2\ijlMR - ^r)4MR + h^)\p, 1; -p, 2) 



8(d) 



— ^(£;) exp(zpr) 



(A12) 



The term proportional to r ^ term in the expansion of Eq. (A12) in powers of r can be matched by the matrix element 
of $1'i>(il) in Eq. (AlO) if the Wilson coefficient of this operator is l/(16'7r^r^), in accord with Eq. (53). 



e. Momentum Integrals 

We list here the momentum integrals that arise in the calculation of matrix elements in the preceding subsections. 
The calculations of the matrix elements of bilocal operators in Sections A 2 a and A 2 d require the following Fourier 
transforms: 

(fq exp{iq-r) 1 /■ ^ /aiq \ 

— ■exp(ipr), (A13a) 



/ 



(27r)^ q^ — p^ — it Airr 

^"P^^*^-") - -^eMipr). (A13b) 



(27r)3 (g2 -p2_ jg)2 g^p 

The calculations of the matrix elements of local operators in Sections A2b and A 2c require the following integrals: 



/ 



+ £, (A14a) 



(27r)3 g2 _ p2 _ 27r2 47r 

d^q 1 i 



{2TTf{q^-p^-ief 8ttp' 



(A14b) 
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The integral in Eq. (A14a) is ultraviolet divergent. It has been evaluated using an ultraviolet cutoff |g| < A. 
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